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Abstract. - We present the quantitative interpretation for the glass-like behavior of thermal 
conductivities k(T) for type-I clathrate compounds involving off-centered guest ions. It is shown 
that the dipole-dipole interaction generated in cage/guest-ion systems is crucial to reproduce the 
characteristics of thermal conductivities for these symmetry-broken clathrates. The above scenario 
also explains well the difference of k(T) between the p-type and the n-type /?-BGS found recently 
by K. Suekuni et al. [Phys. Rev. B, 77 (2008) 235119.] 
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Introduction. — The "phonon-glass electron-crystal" 
concept, proposed over a decade ago [1] in order to ex- 
plore thermoelectric materials with the high figure of mer- 
its Z , has revived interest in thermal transports in solids. 
Certain classes of thermoelectric materials with low ther- 
mal conductivities have been searched, leading to the dis- 
covery of typc-I clathrate compounds with the formula 
IIgIIIi6lV3o [2], where the group III and IV elements con- 
stitute framework atoms of cages and the group II ele- 
ments are guest ions in the cages. Subsequent work [3-11] 
has experimentally confirmed that the thermal conductiv- 
ities k(T) of these compounds are not only drastically sup- 
pressed but also display characteristics identical to those 
of network glasses [12]. 

Avila et al. [11] have achieved the lowest recorded k(T) 
among clathrate compounds in /3-Ba§Gai6Sn3o (/3-BGS). 
In addition, Suekuni et al. [10] have recently discovered 
the clear difference of k(T) between the p-type and the 
n-type of /3-BGS. Actually, they have found that K n (T) 
of the n-type is smaller than that for k p (T) of the p-type 
below a few K, while the reverse tendency in the regime 
above a few K. 



( a ) E-mail: knshtSyukawa. kyoto-u. ac . jp 
^The dimcnsionlcss parameter ZT is often used to evaluate the 
thermoelectric efficiency at temperature T. Here Z = S 2 /p(K p + K e ) 
with p the electrical resistivity, re e (k p ) the electric (phonon) ther- 
mal conductivity, and S the thermopower or Seebeck coefficient, 
respectively. The low n p is crucial to increase the efficiency of ther- 
moelectric conversion since n e is combined with p. 



It is noteworthy that thermal conductivities almost 
identical to those of glasses are observed in the clathrates 
with off-centered guest ions (symmetry-broken clathrates) 
but not in those without off-centered guest ions (symmet- 
ric clathrates) with translational invariance [2-11]. We 
illustrate a scheme of the characteristics of n(T) in Fig. 1 
for both of the clathrates with and without off-centered 
guest ions. Figure 1 shows that n(T) varies as T a (a w 2) 
at temperatures below a few K. Above a few K, there ap- 
pears a plateau in k(T), and it lasts up to 10 K. Above 
10 K, k(T) rises again T-linearly. From 100 K, k(T) curls 
over. These are almost identical to those of general glasses 
with topological disorder [13], while k(T) of the symmet- 
ric clathrates, depicted by a dotted line in Fig. 1, shows 
identical thermal conductivities to those of crystals with 
translational symmetry. 

In this connection, we should mention that orien- 
tationally disordered crystallines, so-called orientational 
glasses, show all of the universal properties of structural 
glasses [14,15]. Grannan et al. [16] and Randeria and 
Scthna [17] have theoretically investigated the role of long- 
range interaction between elastic dipolcs randomly dis- 
tributed in space. Although their systems are related with 
the present work, electric dipolcs in our systems are reg- 
ularly distributed and the situation is clearer. Recently, 
Parshin et al. [18] have discussed the vibrational instabil- 
ity of weakly interacting harmonic modes controlled by the 
4th-order anharmonicity of the potential function. They 



E. Kaneshita et al. 



have demonstrated that the instability produces both of 
features observed in glasses: the so-called boson peak and 
two-level tunneling systems. Our work takes into account 
the explicit interaction between electric dipoles in nearby 
cages for symmetry-broken clathrate systems rather than 
the harmonic coupling between the oscillators. 

In this Letter, we formulate the thermal conductiv- 
ity of the symmetry-broken clathrates in terms of the 
intcracting-dipole picture and provide a theoretical inter- 
pretation on why the symmetry-broken clathrates show 
the characteristics of k(T) almost identical to those of 
network glasses, by focusing our attention on the re- 
cent experiments for the prototype type-I clathrate, (3- 
BGS [8-11]. Our concepts are general and not only appli- 
cable to /3-BGS but also to other types of clathrates with 
off-centered guest ions. 




(i) l (ii) l (iii) l (iv) 
Temperature 

Fig. 1: (Color online) Scheme of thermal conductivities n(T) 
of type-I clathrates. The solid line shows the glass-like temper- 
ature dependence of fc(T) for symmetry-broken clathrate with 
off-centered guest ions. The dotted line gives the crystallinelike 
behaviors of symmetric type-I clathrate without off-centered 
guest ions. In the low-temperature region T < 1 K (i), the 
former behaves as T 2 , while the latter does as T 3 . The plateau 
temperature region appears in the former case at ~ 10 K. The 
thin dashed line depicts the T-linear behavior in the tempera- 
ture region above 100 K (iii), which approaches zero as T — > 0. 
From 100 K, «(T) curls over. 



Thermal conductivities below a few K. The 

positional symmetry of the guest ions in type-I clathrates 
is broken by increasing the size of the cages, and the guest 
ions take off-center positions at ro apart from the center of 
the cage [3-6,8-11]. A typical distortion r = 0.43 A in /?- 
BGS has been determined from diffraction experiments [8, 
11]. This value is about 7.2% of the distance between the 
neighboring 14-hedron cages (a/2 = 5.84 A, where a is the 
lattice constant). 

The deviation of a guest ion from the center of the cage 



in /3-BGS induces a large electric-dipole moment, which is 
due to the difference of the charge between the Ba 2+ guest 
ion and the Ga~ ions of the cage (see Fig. 2) as pointed 
out in [19]. The strength of the electric dipole moment 
becomes p ss 4.1 Debyes for /3-BGS. Thus, the effect of 
the dipolar interaction between guest ions becomes crucial 
since the coupling energy of the nearest neighbor dipole- 
dipole interaction is of the order of 60 K [19]. 



(a) 



(b) 




Fig. 2: (Color online) (a) A 14-hedron cage consisting of anions 
and a guest cation, (b) The cage (outer circle) and the off- 
centered guest ion compose an effective electric dipole moment 
(thick arrow), which is the vector sum of each dipole (thin 
arrow) . 



Of importance is that the long-range and anisotropic 
dipole-dipole interaction constitutes many potential min- 
ima in the configuration space of multiple dimensions; 
therefore, the quantum states should be represented in 
a configuration space. There should be many stable con- 
figurations, and we represent the i-ih stable configuration 
by Xi in the configuration space of multiple dimensions. 
The stable states should execute the zero-point oscillations 
(quasi-ground states) around the stable configurations. 

When two stable states |Xj,{0}) and \Xj,{0}) exist 
near each other in the configuration space, the two-level 
tunneling states \(Xi, Xj), ±) are constructed from these 
states, where + (— ) denotes the higher (lower) energy 
state. The conceptional illustration of the tunneling states 
in the configuration space is shown in Fig. 3 (a). The point 
is that the tunnclings occur in the configuration space of 
multiple dimensions but not in a single cage. 

The thermal conductivity is expressed as 
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where tibe is the Bose-Einstein distribution function, and 
v e , T e , and D e (w) are the group velocity dio{k) / dk, life- 
time, and density of states of thermal phonons of the mode 
e, respectively. We take 
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where 



Fig. 3: (Color online) Scheme of states in configuration space. 
The states for different dipole-configurations Xi and Xj are 
expressed by G (= |X;,{0})) and G (= \Xj,{0})). The hy- 
bridization of these states creates two-level tunneling states T± 
(= \(Xi,Xj), ±)). (b) The excited states \Xi, {n^}) (labeled as 
Af , Bf , etc.) at a certain configuration Xi is illustrated. The 
transition between the excited states Af and Bf corresponds to 
a hopping of the quasilocalized oscillation in the real space (See 
the text). 



where uj e is the upper limit of the extended acoustic- 
phonons with linear dispersion. For v € , we take the average 
v s = 2400 ms" 1 obtained from the observed velocities for 
/3-BGS [10]: v\ = 3369 ms -1 for the longitudinal phonons 
and vt = 1907 ms -1 for the transverse phonons. Note 
that the Matthiesscn's rule l/r e = X^l/ r <M should hold 
when there are different scattering processes. 

In the low-temperature regime below a few K, only two- 
level tunneling states \(Xi,Xj),±) become relevant. The 
transition between {(X^, Xj), — ) and \(Xi, Xj),+) is in- 
duced by a resonant absorption or emission of the thermal 
phonons, and this interaction yields the scattering rate of 
the form 
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where g is the deformation coupling constant, and p the 
mass density (p = 6010 kgm -3 for /3-BGS). e r is the di- 
electric constant, and rj equals to ^tun/^di, where ritun and 
ndi are the numbers of the tunneling states and the actual 
dipoles per volume, respectively; i.e., 1/rj is equivalent to 
the averaged number of off-centered guest ions involved in 
a dipole rearrangement caused by tunneling [19] and is a 
function of p. It should be noted that the analysis of the 
specific heats for /3-BGS provides the average number of 
off-centered guest ions 1/rj ~ 20 [19], implying the failure 
of independent guest ion picture. 

By combining the above scattering rate with Eq. (1), 
k(T) of the type-I clathrates becomes 
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The convenient form with the specific heat Ctun is ob- 
tained as 
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From the experimental data of k(T) ~ 0.02T 2 
Wr 1 nr 1 and C(T) ~ 30T raJmol" 1 FT 1 below a few 
K for symmetry-broken /3-BGS [10], the deformation cou- 
pling g is estimated as g ~ 0.2 eV {g/ks — 2100 K). This 
is a reasonable value because the deformation coupling 
constants in glasses are in the range of 0.1-1.0 eV [20,21]. 

It is worthwhile to consider the effect of the strength 
change of electric dipole moment p on k(T). Suckuni et 
al. [10] have observed the difference of thermal conduc- 
tivities between the n— type and the p— type /3-BGS. The 
observed thermal conductivities below a few K of the n- 
type /3-BGS are a few percent smaller than those of the 
p-type [10]. This difference of k(T) should arise from the 
difference of the dipole moments: The deviation of the 
guest ions is different between the n-type (ro = 0.434) 
and the p-type (0.439 A). This effect appears in Eq. (4) as 
the factor rj/p 2 . (Note that -q should also depend on p.) 

Since the dipolar interaction is proportional to p 2 /r 3 , 
where r is the distance between the dipoles, the effective 
radius r e of the dipolar interaction is proportional to p$ . 
Therefore, the average number of dipoles involved in the 
tunneling, 1/rj, should be proportional to p 2 , then it fol- 
lows that p 2 /rj oc p A or 
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Thus, the difference of the dipoles, p p > p n , leads to the 
relation k p > K n below a few K, while the observed specific 
heats behave the other way, C tU n,p(T) < C tU n,n{T) [10]. 
Both characteristics arc consistent with our formulation 
in Eq. (5), where C tun appears in the denominator. Our 
formulation reproduces the different tendency between the 
specific heats and the thermal conductivities for the n-type 
and the p-type /3-BGS. 

It should be noted that the above formulations are 
valid only when the system can be well-described by the 
interacting-dipole picture. When p is very small, Eqs. (7) 
and (3)-(5) are not valid. In the above formulation, we as- 
sumed the energies of the tunneling states are uniformly 
distributed, which should be true only for the states dis- 
tributed in a wide energy range. For a small p, however, 
thv staters are distributed in a narrow energy range or do 
not exist. 

For such a small p, the scattering rate by the tunnel- 
ing states must become small as the dipolar feature weak- 
ens, i.e., r t ~n — > as p — > 0. Then, the total scattering 
rate of thermal phonons is given by r _1 = r t ~ n + t^ u , 
where (= const.) is the boundary scattering of the 
system and constant. This yields a constant r _1 , leading 
the same Independence of n(T) at low temperatures as 
observed in symmetric clathrates. Thus, k(T) becomes 
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more crystalline-like as p decreases. The temperature de- 
pendence changes from T 2 to T 3 when p gets close to zero, 
as observed experimentally [10]. 

Plateau of k(T) above a few K. — Observed ther- 
mal conductivities n{T) for /3-BGS [10] show the plateau 
behavior above ~ 3 K and increase again T-linearly above 
~ 10 K. In addition, the experiments [10] show the reverse 
tendency K p < K n at the plateau, compared with the case 
of k p > K n below ~ 3 K. The Rayleigh scattering is negli- 
gible owing to the irrelevance of mass-density fluctuation 
in the length scale of excited-phonon wavelengths below a 
few K. 

These observations indicate that other types of oscil- 
lating states become relevant: for example, angular os- 
cillation modes of dipoles as in the case of oricntational 
glasses [14-17]. The eigenfrequency of angular oscilla- 
tion mode, cog, should exhibit the moment-of-inertia ef- 
fect u>e cx I~ x / 2 with the definition / = m^rg, whose 
effect has been clearly demonstrated for oricntational 
glasses [22]. For /3-BGS, the moments of inertia pro- 
vide the inequality I p > I ni which leads the inequality 
ujg^p < (jQe <n for p-type and n-type /3-BGS. The coupling of 
acoustic phonons with angular oscillation modes composes 
the flattened and anti-crossing dispersion curve. 

The expression of k(T) in Eq. (1) suggests that, if the 
scattering rate of thermal phonon r _1 is independent of 
T, k(T) takes the same temperature dependence as the 
Dulong-Petit limit of the Debye specific heat in the high- 
temperature region, i.e., k(T) = const. As demonstrated 
in Rcf. [23], the dominant phonons at fcgT > tiu/ e , namely, 
the coupled modes behave as strongly localized modes sat- 
isfying the condition uigT{ujg) ps 1. Under the situation 
ksT > TiCug, k{T) behaves oc uig from Eq. (1). Hence, 
we obtain the plateau feature in addition to the type- 
dependence, k p (T) < K n (T). 

In this connection, we should mention the experimen- 
tally observed broad peak in specific heat C(T)/T 3 at 
around 5 K for /3-BGS [8—11]. Similar broad peaks are 
observed in all glasses (for network glasses, see, e.g., 
Rcf. [13]) and are called the boson peak. The boson peak 
in the specific heat refers to an excess contribution of 
an excitation over the usual Debye density of states [13]. 
It is natural to expect that the boson-peak-like spectra 
for symmetry-broken clathrates should also be observed 
in spectroscopic measurements. Such peaks have been 
surely found in Raman scattering measurements for /3- 
XGG (X=Ba, Sr, Eu) [24,25]. /3-EGG, a Symmetry- 
broken clathrate, shows boson-peak-like modes in the Ra- 
man scattering spectra with a broad peak around a few 
10 cm" 1 [24,25]. 

In the case of silica glass(u-Si02), hyper- Raman scatter- 
ing experiments performed by Helen ct al. [26] have defi- 
nitely shown the existence of hyper-Raman active modes 
at slightly lower energies than those of Raman active 
modes. Thus, there should exist corresponding additional 
modes (angular-oscillation modes) with broad spectra for 



symmetry-broken clathrates. These are weakly overlapped 
with Raman active modes as in the case of w-Si02- The 
existence of these modes should be confirmed by hyper- 
Raman, infrared, or inelastic neutron scattering experi- 
ments. 

We should mention that the symmetric clathrates, i.e., 
the type- VIII clathrates or type-I clathrates without off- 
centered guest ions (e. g. a-BGS; /3-SrgGai6Si3o, etc.) 
show similar peaks in the specific heats at higher temper- 
atures [9,10]. In these symmetric compounds, the excess 
specific heats have been interpreted in terms of a har- 
monic oscillator, i.e., the Einstein model [9,10], unlike for 
the symmetry-broken clathrates. 

T-linear rise above plateau region. At temper- 
atures above the plateau region, the thermal conductivi- 
ties follow the relation n(T) = 7T [3-11]. Actually, the 
factor 7 takes 7 = 0.009 WK _2 m _1 for the n-type and 
7 = 0.007 WK^m- 1 for the p-type /3-BGS [8-11]. Of 
importance is that these behave as n(T) — > when extrap- 
olating the temperature as T — > (see Fig. 1), although 
this feature has not been pointed out so far. It implies that 
a new additional heat-transport channel opens up above 
the plateau region [27,28]. The additional channel is as- 
sociated with the hopping of the quailocalizcd oscillations 
kicked by the thermal phonons as shown below. 



A 




A' 




Fig. 4: The hopping processes of quasilocalized modes. A ther- 
mal phonon (wavy lines) kicks a quasilocalized mode (solid 
lines) . 



We consider an interaction where a thermal phonon 
with wave vector k kicks a locally oscillating state, 
\Xi,{n\}), into another, \Xj,{n' x }), which corresponds 
to the A1<->B1 transition in Fig. 3 (b). The interaction 
Hamiltonian is expressed by 
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where the operators b^ and c 1 (b and c) are the creation 
(annihilation) operators of the thermal phonons and the 
quasilocalized modes, respectively. The coefficient r e ff is 
the coupling constant of the three phonon processes shown 
in Fig. 4. According to Ref. [28], this type of process yields 
a T-linear contribution to the thermal conductivity. Using 
the definition «hop (T) = Ea #aMo P K, T) and following 
a discussion similar to that in Ref. [28], the hopping con- 
tribution to the thermal transport turns out to be 



Khop(T) = 



i2 2 r c 2 ff fc|r 
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where p is the mass density, which is 6.01 x 10 3 kgm -3 for 
/3-BGS. Here, we defined (4ir/3)£ 3 as the average volume 
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for finding a single quasilocalizcd mode. Since the con- 
tribution of the thermal phonon to the total n is reduced 
owing to the scattering process by the hopping interaction 
between the low-energy thermal phonon and the quasilo- 
calizcd modes (local vibrations bilinearly coupled to the 
acoustic phonons), wc assume k(T) ~ Khop(T) and obtain 
the T-linear thermal conductivity. 

The localization length £ is determined so that 
(4ir/3)£ 3 rj (l/6?7)a 3 , where 6 is the number of dipoles in 
the unit cell. Here, we set r\ v ~ 0.060 and r) n ~ 0.58ry p for 
/3-BGS from Eq. (8) and the speculation in Ref. [19]. From 
Eq. (10) and the experimentally observed k(T) (oc 7T), 
the coupling constant T e g is estimated as \T e g | ~ 3.1 x 10 13 
(2.7 x 10 12 ) for the n-type (p-type) /3-BGS. The estimated 
coupling constants are of the same order as in the case 
of network glasses [29], where the large anharmonicity is 
the necessary condition for explaining the T-linear rise of 
k(T) above the plateau [28]. 

We mention the numerical simulations of thermal con- 
ductivities for amorphous Si [30] performed under the fol- 
lowing conditions: The system size was 10 3 atoms with 
periodic boundary conditions and a harmonic approxima- 
tion. The finite size of 16.3 A prohibits any lower-energy 
propagating modes below 10 meV, which indicates their 
results are applicable at higher temperature regime ob- 
served for amorphous Si by Cahill et al [31]. However, 
their numerical results do not produce the extrapolating 
behavior, re(T) — > at T — > 0, experimentally established 
for glasses and symmetry-broken clathrates. 

Propagating acoustic modes at very low-energies do ex- 
ist in glasses at temperatures and play a role for thermal 
conductivity above the end of the plateau since they obey 
the Bose-Einstein distribution. Our thermal-conductivity 
mechanism above the plateau highlights the anharmonic 
interaction between quasilocalizcd modes and these prop- 
agating acoustic phonons incorporating the strong anhar- 
monicity in symmetry-broken clathrates compared with 
the case of amorphous Si. 

Summary. — We have provided the interpretation 
of the glass-like behavior of the thermal conductivities 
of symmetry-broken clathrates over a wide temperature 
range. Considering the scattering of the thermal phonons 
by the tunneling sates and the states coupled with angu- 
lar orientation modes, in addition to the hopping of these 
quasilocalized modes, we have reproduced the tempera- 
ture dependence of k(T) and have obtained the following 
conclusions. 

i) Two-level tunneling states generated in a configura- 
tion space under the interacting-dipole picture explain well 
the observed T 2 -dependence of k(T) below a few K, and, 
in addition, the feature k p > n n between the p-type and 
n-type symmetry-broken clathrates below a few K. 

ii) The plateau of k(T) is due to the coupling of acoustic 
phonons with angular orientation modes of dipoles [16], 
namely, the crossover from extended to quasilocalized 
phonons (the flattening of dispersion curves of thermal 



phonons) should occur in the acoustic phonon branch. The 
observations n p < n n [10] at the plateau region arise from 
the difference of angular oscillation, ug >p < u)g tn . We have 
predicted the existence of low-energy modes at around 
10 cm -1 for /3-BGS attributing to the angular oscillation 
modes of dipoles, which are active to infrared or hyper- 
Raman scattering. 

iii) The observed k(T) above - 10 K for /3-BGS be- 
haves as k(T) — > as T — > 0. This suggests the onset of a 
new additional heat transport channel. We have proposed 
the new channel associated with the hopping of quasilo- 
calizcd modes kicked by the thermal phonons, namely, the 
hopping mechanism of the quasilocalized modes comes up 
above the plateau temperature region and leads to a heat 
transport with the T-linear dependence of k(T). 
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